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ENERGY EXCHANGES, RATE COEFFICIENTS AND SPECTRAL DISTORTIONS 
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Energy exchange processes play a crucial role in the early Universe, affecting the thermal balance and 
the dynamical evolution of the primordial gas. In the present work we focus on the consequences of a 
non-thermal distribution of the level populations of H2: first, we determine the excitation temperatures 
of vibrational transitions and the non-equilibrium heat transfer; second, we compare the modifications 
to chemical reaction rate coefficients with respect to the values obtained assuming local thermody- 
namic equilibrium; third, we compute the spectral distortions to the cosmic background radiation 
generated by the formation of H2 in vibrationally excited levels. We conclude that non-equilibrium 
processes cannot be ignored in cosmological simulations of the evolution of baryons, although their 
observational signatures remain below current limits of detection. New fits to the equilibrium and 
non-equilibrium heat transfer functions are provided. 

Subject headings: molecular processes; cosmology: early Universe, cosmic microwave background 



1. INTRODUCTION 

Understanding the thermal evolution of the Universe 
in the epoch in which atoms and molecules formed is a 
crucial step to properly model the birth of the first bound 
structures (e.g.. Flower & Pineau des Forets 2001). In 
particular, the balance between cooling and heating pro- 
cesses has to be taken into account and modeled accord- 
ing to the chemical and physical processes occuring in the 
primordial plasma. It is well established that Lyman- 
a cooling is effective at gas temperature higher than 
~ 8000 K, corresponding to redshifts z > 2700, while 
primordial molecules, in particular H2 and HD formed 
at 2; < 1000, are the most efficient cooling agents of the 
pristine plasma at lower temperatures. 

Several authors have calculated the heating and 
cooling fun c tions of the primordial molecul ar species: 
iPalla et all (Il98l. Lep p fc ShuU (1984), IP^ et al.l 
([1991 . iLe Bourlot et all ([T999h. IP^fc Signord ([1991 
Galli & Palla 1998 (hereafter GP98'). iCoppola et al.l 
(|2011bf ). One of the standard assumptions in these cal- 
culations is that the population of internal states can be 
described by a Boltzmann distribution. This hypothesis 
is valid in many astrophysical environments where the 
density is sufficiently high to bring the internal degrees 
of freedom to a condition of local thermodynamic equilib- 
rium (LTE). For a given species, e.g. H2, this condition is 
quantified in terms of a critical density ncr(H) defined as 
the ratio between radiative and collisional de-excitation 
coefficients of H2 . It is easy to check that even at highest 
redshifts (z « 1000), the ambient baryon density rib is 
~ 2 orders of magnitude below the critical density ricr (H) 
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(although considering specific rotational transitions the 
critical density is below the baryon one up to lower red- 
shift z « 500). In addition, several physical phenomena 
(e.g. shocks) and chemical processes can produce de- 
viations from LTE. In particular, most of the gas-phase 
molecular formation processes selectively produce species 

in states that deviate significantly from L TE. 

In th e case of the early Universe, ICoppola et al.l 
(|2011a( ) (hereafter Cll) computed the vibrational distri- 
bution of H2 and HJ formed at redshifts 10 < z < 1000 
and found that high sovrathermal tails are present es- 
pecially at low z. The existence of non-equilibrium fea- 
tures in the level populations of H2 is important, be- 
cause of the role of this species as a coolant of primor- 
dial gas. In the present work, we extend the work of Cll 
to study the following physical quantities relevant to the 
non-equilibrium energy exchange in the primordial Uni- 
verse: excitation temperatures, heat transfer functions, 
reaction rates and spectral distortions of cosmic back- 
ground radiation (CBR). 

2. REDUCED MODEL: ORTHO- AND PARA- STATES 

Because the number of roto- vibrational levels involved 
in the kinetics of H2 is too high (w 300) for a direct ex- 
tension of the approach used by Cll, in the present work 
we implement a reduced model in order to provide a sim- 
pler starting point for more extended calculations. This 
model is based on the assumptions that: (1) the most 
important channels determining the population of vibra- 
tional states are the associative detachment reaction 

H- + H ^ H2 + e, 

and the spontaneous and stimulated radiative transi- 
tions between roto-vibrational levels; and (2) a steady- 
state approximation can be applied to the kinetics of 
vibrational levels. Under these hypotheses, two steady- 
state Master Equations (for ortho- and para- states, re- 
spectively) can be written for the 15 vibrational levels 
{i = 0,14) of H2: 

ft ^ Rij - ^ Rj 
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where fi is the fractional abundance of H2 in the 
level, /h- and /h those of H~ and H respectively, Rij 
is the matrix of radiative coefficients including absorp- 
tion processes, calculated as in Cll averaging over the 
initial rotational levels and summing over the final ones 
the E instein coefficients computed by IWolniewicz et all 
(|1998D . The values of /h- and /h are taken from the 
complete kinetic model by Cll. The rate coefficients fc,; 
for the associative detachment reaction for the i^^ vi- 
brational level formation have been evaluated using the 
cross-sections CTy calculated bv ICizek et al.l ()1998[ ) sum- 
ming over all final rotational states: 



i,„ax(») ^ 

ki (Tin ) — - ^ ^ (Tin ) , 

3=0 



(2) 



where 



Tr^i{k^TmY Jo 



E<7,j{E)e 



-E/{k-BT,„ 



dE, 



(3) 

with /i reduced mass of the system and T„i matter tcm- 
pe rature. It should be noted that experimental results 
bv lKreckel et al.l ()201Cl[ ) recently showed very good agree- 
ment with these quantum calculations. 

The equation for the vibrational ground state /o in 
Eq. [T] is replaced by the normalization condition 



fi — /h2 , 



(4) 



where /h2(-2^) is the fraction of H2 which is also taken 
from the complete model. Figure [1] shows the fractional 
abundance of vibrational levels obtained using the re- 
duced model described in the present section. The results 
for all fi are satisfactory close to those obtained with the 
fully kinetic model shown in Figure 10 of Cll, at least 
for not too high values of z. These results show that the 
hypothesis of steady-state can be applied to the present 
problem with enough confidence to proceed to the study 
of the ortho- and para- states. For this we calculate two 
sets of rate coefficients fci_ortho(Tni) and fc,;.para(Tm) for 
each vibrational level i, 



^2,ortho(Tni) ^ ^ ^zj(Tjn), 



■j odd 



(5) 



The Rij coefficients are also thermally averaged over a 
partial distribution. The first equation of each system is 
replaced by the normalization condition 



itho/pa 



1, 



(6) 



which means that we are calculating the vibrational dis- 
tribution of each of the two species but not the total 
fraction of ortho- and para- hydrogen, which cannot be 
calculated by a steady-state approach. 

3. EXCITATION TEMPERATURES 

For each transition 0-w, the excitation temperature is 
defined as 

_ E^ - Eq 
J-O-v — 7—, -, T, [I) 

kB in{no/ny) 




Fig. 1. — Fractional abundances of the vibrational levels of H2 
according to the reduced steady-state model described in the text. 

where Ey and riy represent the energy and fractional 
abundance of the w*^ vibrational level, respectively, and 
fee is the Boltzmann constant. 

In Figure [21 the excitation temperature of the transi- 
tions — 1) is compared to the temperatures of matter 
and radiation as a function of z. For better clarity, a few 
curves have been repeated in the two panels. The fig- 
ure shows that vibration decouples from the translational 
motion, and this happens at higher z than the decoupling 
of radiation and matter; for the most excited levels this 
corresponds to an epoch where the conditions are near 
equilibrium. Later, after the decoupling of matter and 
radiation, To_i is slightly larger than the radiation tem- 
perature Tj-; the excitation temperatures Tq-u for higher 
V are progressively higher. This result can be explained 
by considering the reaction network: the vibrational lev- 
els are formed by associative reactions that preferentially 
populate highly excited levels, while the vibrational man- 
ifold as a whole is coupled to the CBR (which provides 
a heat sink) much better than to the matter; this lat- 
ter coupling occurs via the relatively ineffective II2/II 
VT processes. Therefore, all levels pairs are expected to 
be warmer than the radiation field, but the lowest pairs 
are closer to equilibrium with the radiation because the 
chemical heating is lower. 

The separation of To_i and occurs at z « 100, which 
brings this phenomenon not far from potential indirect 
observation (e.g. effects on reaction rates of processes 
and consequent different fractional abundances of chem- 
ical species at lower z). Another relevant feature of our 
calculation is that To_i is stable at about 200 K at the 
age of formation of the first structures. Since the 1 — >■ 
transition is an important heat radiator, our results in- 
dicate that To_i is a more appropriate initial condition 
for the vibrational temperature of H2 in hydrodynamic 
collapse models than T^, which is considerably lower. In 
Table [T] the values for z^ec (redshift at which excitation 
temperature decouples from radiation temperature) and 
Zfreeze-out (rcdshift at which the freeze-out temperature 
is reached) are reported for each To-i] the former are 
evaluated considering a relative deviation from Tj- larger 
than 1%, while the latter are calculated searching for rel- 
ative deviation of To-i at each z from the value at z = 10 
smaller than 10^"^. 

Figure [3] shows the values of the relative difference be- 
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Fig. 2. — H2 excitation temperatures Tq-u, compared to the tem- 
perature of the radiation {solid curve) and matter {dashed curve). 
Top panel: d = 1 to d = 5; bottom panel: v = 6 to v = 14. 



TABLE 1 

Decoupling redshift of excitation temperatures 



Excitation temperature (K) z^^^ 



rcczc — out 



7b- 14 


697 


382 


7b- 13 


689 


379 


To- 12 


676 


377 


7b-ii 


647 


379 


7b- 10 


627 


382 


7b-9 


593 


537 


7b-s 


551 


524 


7b_7 


511 


389 


7b-6 


457 


389 


To-5 


391 


377 


To-4 


363 


346 


To-3 


300 


313 


To-2 


198 


174 


To-1 


108 


84 



tween the vibrational excitation temperatures of ortho- 
and para- states. As can be seen, significant deviations 
between vibrational temperature of states with different 
rotational symmetry can be detected at low z and high 
i. These differences suggest that the issue of rotational 
non-equilibrium could be important and deserves to be 
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Fig. 3. — Ratio of excitation temperatures for ortho- and para- 
states given as (To-i^ortho - To-i,para)/To-i, ortho ■ Top panel: v = 
1 to 1; = 8; bottom panel: v = 9 to v = 14. 

addressed accurately in future studies, i.e. by solving the 
Master Equation for a full rotovibrational manifold. 

4. HEAT TRANSFER FUNCTION 

In this paper we discuss the role of chemical energy 
from exothermic reaction, which is ultimately dissipated 
either into radiation or into the thermal energy of H 
atoms. This energy flow is described by heating and 
cooling functions, usually indicated with the symbols F 
and A, respectively. We consider the net molecular heat 
transfer, defined as the sum of all radiative excitations 
of H2 followed by collisional de-excitations with H atoms 
and all collisional excitations followed by radiative decay: 



$(rn„Tr) = (r-A)H, = 



1 



E 

(v' ,j')<{v,j) 



(8) 



n, 



where Tm and Tj- are the temperatures of matter 
and radiation, k(^^i jt)_:f(^y j--j is the VT (vibrational- 
translational) rate coefficient for H-|-H2(w', j') — > 
H+H2(w,j); and describes the distribution of ro- 
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tovibrational levels, 



gjnv{2j + l)exp I - 



-E^.■J--E^■,0 



(9) 



with Qj equal to 1/4 and 3/4 for the para- and ortho- 
states, respectively, and Zy{Tr) rotational partition func- 
tion for the w*'' level. It can be seen from Eqs. ([S])-® 
that the rotational energy is distributed according to 
the Boltzmann law with temperature Tj. while VT co- 
efficients depend on the matter temperature Tm. Thus, 
the heat transfer function depends on two temperatures, 
as a consequence of the different couplings of the internal 
degrees of freedom of molecules with the radiation and 
the matter: a faster coupling occuring between radiation 
and rotation, a slower coupling between translation and 
vibration. 

We explore both equilibrium and non-equilibrium 
cases, corresponding to different values of n.„; in the for- 
mer case the vibrational levels are distributed following 
the Boltzmann population equation, whereas in the latter 
we adopt the level populations resulting from the kinetic 
model of Cll. VT rate coe fficients have been taken from 
lEsposito etall (Il999l[200l . These coefficients arc given 
as functions of the initial and final rotovibrational quan- 
tum numbers and of temperature allowing the inclusion 
of the full sets of collisional transitions in our calcula- 
tion. The equilibrium one-temperature heat transfer us- 
ing the same rate coefficients has been compared with 
the analytical expression for the cooling function given 
by GP98. In GP98, the gas was not embedded in any ra- 
diation field; consequently, only collisional dc-excitations 
were considered. 

As it can be seen from Figure HI the non-equilibrium 
heat transfer function decreases more rapidly than the 
equilibrium one at lower temperature, because of the in- 
creased efficiency in the energy exchange due to the long 
sovrathermal tails in the vibrational distribution. In the 
two temperatures calculation, the deviation is amplified 
as a consequence of the strong decoupling in the energy 
exchange among degrees of freedom and of the differ- 
ent trend of the gas and radiation temperatures. The 
first effect, which is due to the deviation of the vibra- 
tional population from the equilibrium distribution (es- 
sentially those of the lowest levels) is seen at z « 1000, 
and amounts to a factor of ~ 2. The effect of the sep- 
aration of Tm and occurs at low z where these two 
temperatures are considerably different. Consequently, 
in the computation of the heat transfer function at least 
two temperatures must be used: one is the temperature 
of the level population and the other is the translational 
temperature. This difference is evident in Eq. 

All heat transfer functions available in the literature 
are calculated assuming a single temperature, usually 
set equal to T^. Such usage cannot capture the second 
non-equilibrium effect described above, since ^^('T^ni, X^-) 
is implicitly set equal to k{Ti). A better solution is to 
use in the context of early Universe models heat transfer 
functions calculated including non-equilibrium effects al- 
though fitted later as a function of a single temperature. 
Fits for the non-equilibrium case using the present two- 
temperatures model and the equilibrium case using the 
newest available data by lEsposito et al.l ()1999l |2001[ ) are 
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Fig. 4. — Heat transfer function (F — A)h2 ^ ^ function of the 
radiation temperature Tr. Dotted line: GP98 fit for the cooling 
function; solid line: LTE calculation with present VT coefficients; 
dashed line: non-equilibrium vibrational distribution contribution. 
Calculation are reported both in the one- temperature case {top 
panel) and the two-temperatures one {bottom panel). 



TABLE 2 

H2 HEAT TRANSFER FUNCTION 



Fitting coefficients 


equilibrium 


ao 


= -145.05 




ai 


= 136.085 




12 


= -58.6885 




13 


= 11.2688 




a4 


= -0.786142 


non-equilibrium 


ao 


= -393.441 




ai 


= 588.474 




12 


= -380.78 




13 


= 123.858 




04 


= -20.1349 




as 


= 1.30753 



obtained in the form: 



N 



(10) 



n=0 



The coefficients an are listed in Table [21 For the non- 
equilibrium case, the validity of the fit is up to T « 100 K 
(additional data are available upon request). 
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Fig. 5. — Rate coefficients as a function of 2: LTE approxima- 
tion {dashed line, fit by Coppola et al. 2011a) and non-equifibrium 
vibrational distribution function {solid line). 

5. NON-EQUILIBRIUM REACTION RATES 

Using the real non-equilibrium vibrational distribu- 
tions, vibrationally resolved rate coefficients have been 
recomputed and compared with the corresponding LTE 
fits by Cll. In Figure [5] the results for the following pro- 
cesses introduced in the model are shown, both for H2 
andHj: (1) H2(i;)/H+ charge transfer (2) H2(w)/e- dis- 
sociative attachment (3) H^(w) dissociation by collisions 
with H (4) il2{v) dissociation by collisions with H+. The 
LTE fit for dissociative attachment of H2 is taken from 
CapiteUi et al. (2007). Strong deviations from the LTE 
fits can be noted, due to the non-equilibrium pattern, 
especially at low z, where the hypothesis of Boltzmann 
distribution of the vibrational level manifold fails. The 
peak at z « 300 corresponds to that on (and con- 
sequently to H2, via the process of charge transfer with 
H"*"). It should be noted that, in the case of H2 disso- 
ciative attachment, the non-equilibrium calculation fol- 
lows the trend reported by CapiteUi et al. (2007), where 
a simplified model for the non-equilibrium distribution 
was assumed. 
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6. SPECTRAL DISTORTIONS OF THE CBR 

The photons created in the H2 formation process pro- 
duce a distortion of the black-body spectrum of the cos- 
mic background radiation (CBR). Since the maximum 
production of H2 occurs in vibrationally excited states 
by associative detachment in collisions of H and H~ at 
redshifts below z « 100 (see Coppola et al. (2011a)), 
the emission of rovibrational transitions with wavelength 
A « 2 /xm is redshifted today at A « 100-200 /im, in 
the Wien part of the CBR. An early estimate of this 
distortion, based on the rovibrational-resolved associa- 
tive detachment cross sections computed by Bieniek & 
Dalgarno (1979), was made by Khersonskii (1982) while 

Shchekinov & Entel (1984) developed a model for the 
molecular hydrogen distortion due to secondary heating 
processes. We reconsider here the process of vibrational 
emission of primordial H2 molecules with our updated 
chemical network and with a fully kinetic treatment of 
the level populations of H2. 

For each transition from an upper level Vu to a lower 
level vi , with level populations n„ and n/ and degeneracy 
coefficients and g;, the relative perturbation in the 
CBR at the present time is given by 

^"^^ - [^(Zint) - l]T(zint), (11) 

2=0 



where 

'^'(^int) 



giriuizint) 



exp 



hv. 



ul 



is the source function, and 

T{zint) = Z 



91 



1 



ginu{zint) 
guni{zint) 



Hzizint) ' 



(12) 



(13) 



is the redshift-integrated optical depth (sec e.g. Ap- 
pendix A of Bougleux & GaUi 1997). In Eqs. and 
(|13p . 2int is the interaction redshift, at which the ob- 
served frequency v is equal to the redshifted frequency 
of the transition, i.e. 



v{l + Zint) = l^ul, 



(14) 



Aui are the Einstein coefficients, and Hz is the Hubble 
function 

Hz = Ho[nr{l + z)'^ + n„^{l + zf + nkil + z)^ + nA]^/^ 

(15) 

(see Coppola et al. (2011a) for a definition of the cosmo- 
logical constants and their adopted values). 

Figure[6]shows the emission produced by H2 transitions 
with Av = 1, 2, 3 and 4 in the frequency range v = 10- 
1000 cm~^, corresponding to wavelengths A = 10 fim- 
1 mm. To avoid confusion, only the first 4 transitions 
for each Av are shown (i.e., w = l— >0,2— j-l, 
3^2 and 4 ^ 3 for Av = 1, etc.). The figure also 
shows the CBR in the Wien region, and, for compar- 
ison, the spectral features produced by the cosmologi- 
cal recomb ination of H and He (e.g. Chluba fc Sunyae v 
2007.2008. IChluba. Rubifio-Martin and SuiivaevI (|2007t ). 
iRubifio-Martni. Chluba and SunvaevI (|2008[ )). The lat- 
ter are mainly formed by redshifted Ly-a and two-photon 
transitions of H and the corresponding lines from He 
(Chluba & Sunyaev 2010 have recently produced new 
results for this last contribution). 
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Fig. 6. — CBR spectrum at 2=0 together with 
spectral di stortions due to H and He reco mbina - 
tions llChluba. Rubino-Marti'n and SunvaevI 1120071 ). 

IRubiiio-Martin'^Mub^lmT'Sunvae^ II2008I )) and to non- 
equilibrium vibrational molecular transitions for H2. The 
contribution of multiquantum transitions up to A?; = 4 are shown. 

The difficulty of detecting spectral distortions in the 
Wien side of the CBR, in the presence of an infrared 
background (both Galactic and extragalactic) several or- 
ders of magnitude brighter, have been discussed by Wong 
et al. (2006). While a direct detection appears challeng- 
ing (see also Schleicher et al. 2008), we stress that an ex- 
cess of photons over the CBR at wavelengths shorter than 
the peak could represent a significant contribution to sev- 
eral photodestruction processes, as shown by Switzer & 
Hirata (2005) for the photoionization of Li and Hirata 
& Padmanabhan (2006) for the photodetachment of H^. 
Another possibility is fluorescence, i.e. the absorption of 
the short-wavelengths non-thermal photons by atoms or 
molecules followed by re-emission at longer wavelenghts 
in the Rayleigh- Jeans region of the CBR, as suggested 
by Dubrovich & Lipovka (1995). These issues will be 
addressed elsewhere. 

7. CONCLUSIONS 

We have considered several vibrational non- 
equilibrium effects on the chemistry and physics of 
early Universe. Although our present calculations raise 
several questions about the very use of the concept of 
temperature in the early Universe, the most important 
issue concerning thermal transfer and chemical reactivity 
has been addressed. Our calculations show that the 
differences between the excitation, translation and radi- 
ation temperatures can affect the heat transfer functions 
of important species, an effect here demonstrated for 
H2. Our results underline the necessity to fully include 
the consequences of the temperature separations which 
occur at different epochs in the chemical and physical 
evolution of the early Universe. 

Excitation temperatures appear to be higher than 
radiation temperature at low z; this "chemical" pre- 
heating should be considered while modeling the for- 
mation of galaxies, together with virialization heating 
and other physi c al mechanisms usuall y suggested (e.g. 
iMo et al.l ([20(15h . IWang fc Abell (fT996h V We have also 
assessed the hypothesis of steady-state for the vibra- 
tional distribution presenting a reduced kinetic model, 
and calculated the deviations between vibrational tem- 
peratures of ortho- and para- states as a first step towards 
a full non-equilibrium rotovibrational kinetics. Heat 
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transfer functions are calculated for both equilibrium 
and non-equilibrium cases, considering also a novel two- 
temperatures approach that takes into account the dif- 
ferent rates of energy exchange among molecular degrees 
of freedom. For pure vibrational transitions, the critical 
density is greater than the baryon density, so that the 
hypothesis of non-equilibrium is also valid at higher z. 
Resolving rotations and vibrations gives different results, 
making the limit for z lower (e.g. for the (0, 2) (0, 0) 
transition the critical density is about k, 2.7x10^ m~'^ at 
z « 500.) 

We have evaluated the effects of vibrational non- 
equilibrium on reaction rates. A general increase has 
been pointed out because of the formation of long 
sovrathermal tail in the vibrational distribution, espe- 
cially at low z] this evidence should be added to the 
increase of rate coefficients due to the inclusion of the 
entire vibrational manifold, that by itself can affect in a 
de ep way the fate o f the system modeled (as described 
bv lSethi et al.l (|2010[ ) for the dissociative attachment pro- 
cess) . 



We have computed the spectral deviations to the CBR 
due to the non-equilibrium level populations, considering 
all the transitions. Although the present Planck exper- 
iment and the upcoming James Webb Space Telescope 
(JWST) are able to detect galaxies at high redshift, a 
direct observation of this effect is challenging; for this 
reason, an alternative study of the non-thermal vibra- 
tional photons on the photochemical pathways of atomic 
and molecular kinetic should be undertaken. 
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